CONVERGENCE OF FROZEN GAUSSIAN APPROXIMATION 
FOR HIGH FREQUENCY WAVE PROPAGATION 



JIANFENG LU AND XU YANG 

Abstract. The frozen Gaussian approximation provides a highly efficient 
computational method for high frequency wave propagation. The derivation 
of the method is based on asymptotic analysis. In this paper, for general lin- 
ear strictly hyperbolic system, we establish the rigorous convergence result for 
frozen Gaussian approximation. As a byproduct, higher order frozen Gaussian 
approximation is developed. 



1. Introduction 

This paper is devoted to the proof of convergence of the frozen Gaussian approx- 
imation, introduced in [11,12], for high frequency wave propagation. Numerical 
computation of high frequency wave propagation is an important problem concerned 
in many areas, such as seismic imaging, electromagnetic radiation and scattering, 
and so on. The problem is challenging for direct numerical discretization because 
the mesh size has to be chosen comparable to the wavelength or even smaller in 
order to get accurate solution, however the domain size is usually large so that 
the computational cost is formidably expensive. To look for efficient computation, 
the development of asymptotics-based algorithms has received a great amount of 
attention in recent years. 

The investigations have been focused on two methods: geometric optics and 
Gaussian beam method. The computational methods based on geometric optics 
(see the review articles [3, 19], and references therein) solve eikonal and transport 
equations instead of original hyperbolic system. This makes the choice of mesh 
size frequency-independent, and hence the methods are quite efficient. However, 
eikonal equation can develop singularities that makes the asymptotic approximation 
break down at caustics. To overcome this problem, Popov introduced Gaussian 
beam method in [15], which constructs solution near geometric rays using Taylor 
expansion. Ralston [17] showed that the method gives a valid approximation at 
caustics. One shortcoming of Gaussian beam method is however, since it is based 
on Taylor expansion, the constructed beam solution has to stay near geometric 
rays to maintain accuracy. Therefore the method loses accuracy when the solution 
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spreads [12, 14, 16]. The problem is one of the niajor concerns in application of 
Gaussian beam methods in areas like seismic imaging; see for example [2,5]. 

The frozen Gaussian approximation was proposed in our previous works [11,12] 
to overcome the problems of the aforementioned methods: Geometric optics breaks 
down around caustics; Gaussian beam method loses accuracy when beam spreads. 
The frozen Gaussian approximation is based on asymptotic analysis on phase plane, 
motivated by the Herman-Kluk propagator developed in chemistry literature [4,6, 
7]. It provides a highly efficient computational tool for computing high frequent 
solution to linear hyperbolic system. 

Our previous work [11, 12] developed numerical algorithms for frozen Gaussian 
approximation in both Lagrangian and Eulerian framework. Numerical examples 
indicate the efficiency and accuracy of the method. In the current work, we prove 
the convergence of frozen Gaussian approximation. Denote the propagator of the 
first order strictly hyperbolic system as Vt, and the propagator of the ifth-order 
frozen Gaussian approximation as P/^- The main result of this work is the following 
theorem. 

Theorem. Let {u'^}^^q be a family of asymptotically high frequency initial con- 
ditions, with ||uoIIl2(r3) < M , then we have 

\\Vtul-PlKul\\^,^^,^<e''M. 

Please refer to Section 4 and Theorem 4.1 for the formulation of the frozen 
Gaussian approximation and a precise presentation of the main theorem. The 
asymptotically high frequency initial condition is defined in Definition 2.2. 

Related works. The frozen Gaussian approximation is motivated by the Herman- 
Kluk propagator developed in the chemistry literature [4,6,7], which is used in the 
semiclassical regime of time dependent Schrodinger equation. The convergence 
of Herman-Kluk propagator was recently proved by Swart and Rousse [20] and 
Robert [18]. It is proved that the Herman-Kluk propagator converges to the true 
propagator of the Schrodinger equation. In particular, this means that applied to 
any initial data in L^, the Herman-Kluk propagator provides an accurate result as 
e 0. 

The difference of the Schrodinger equation and the first order hyperbolic system 
lies in the fact that the hyperbolic system might present singularities at p = 
on phase plane (see Section 3 for more details). Therefore, one can not hope to 
obtain similar results for first order hyperbolic system as for Schrodinger equation. 
Indeed, we construct a counter-example which shows that the method fails to give 
a good approximation for low frequency initial data in Example 4.2. On the other 
hand, the results of this work show that frozen Gaussian approximation works for 
the initial data that are high frequent. This is of course the working assumption 
for high frequency wave propagation. 
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The convergence of the Gaussian beam method has been recently investigated in 
[1,8-10]. The results in [1, 10] showed that the K-th order Gaussian beam method 
converges to the true solution with an accuracy of 0{e^/'^). The K-tYi order frozen 
Gaussian approximation, as indicated by Theorem 4.1, has a convergence order 
0{£^). We refer to [12] for a more detailed numerical comparison between frozen 
Gaussian approximation and Gaussian beam method. 

Organization of the paper. In Section 2, we introduce some necessary notations 
and preliminaries for phase plane analysis. The hyperbolic system we considered 
is presented in Section 3. Section 4 describes the formulation of frozen Gaussian 
approximation and states the main convergence theorem. The proof is based on 
the construction of high order approximate solution given in Section 5 and stability 
from the well-posedness of hyperbolic system. We conclude the proof in Section 6. 



2. Preliminaries 

2.1. Notations. We will in general use x, y €E R'' as spatial variables, {q,p) G 
as the variable for the phase space, d is the spatial dimensionality. We will use 
the same notation |-| for absolute value, Euclidean distance, vector norm, (induced) 
matrix norm and (induced) tensor norm. 

For 5 > 0, we define the closed set Kg C M^'' as 

(2.1) Ks = [{q,p) e M^'^l \q\ < 1/(5, \p\ e [6,1/S]}. 
For / : R^d ^ c, we define for fc e N and 6 G R+, 

(2.2) Afc,4/] = max sup 9"/(g,p) 

where a is a multi-index. This definition can be extended to vector valued and 
matrix valued functions straightforwardly. 

For / : ^ and M, N e N, the matrix valued function 9^/ : R^ 
£MxN (Je£ned with the convention that {dxf{x))jk — dx^fkix) for j = 1, • • • , M, 
k = l,---N. 

We use the notations S, C°° and for Schwartz function class, smooth func- 
tions and compact supported smooth functions respectively. 

For convenience, we use the notation 0{£°°): = 0{£°°) means that for any 

lime-'= lA^I = 0. 

e— i-O 

Notations c and C will be used for constants, whose values might change from 
line to line. Sometimes, we will use notations like Ct,k to specify the dependence 
of the constant on the parameters T and K. 
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2.2. Wave packet decomposition. For {q,p) G M. , define ijjf^ ^ as 
(2.3) rg,p{x) = (27r£)-<^/2 g^p|^^p . _ -l\x- /e). 

Define FBI transform on S{W^) as 



(2.4) 



{^^f){q,p) = {^er'''H%,,J) 



= 2-'^/2(7re)-3''/'* / e-'P (^-«^/''-5l^-«l'/7(a;)da;. 

Tlie inverse FBI transform (^^)* on »S(M^'^) is given by 

(2.5) ((^")*5)(a;) =2-'^/2(^g)-3d/4 /" /" g^p^— ^gdp. 

We summarize some elementary properties of the FBI transform, wliose proof is 
standard and can be found, for example, in [13]. 

Proposition 2.1. For any f e S{W^) and g e ^(M^rf), 

(2-6) |l^'^/||i2(]K2d) = ||./|lL2(Rd) ; 

Hence, the domain of and (J^^)* can be extended to L^(R'') and L'^iM?''-) re- 
spectively. Moreover, we have 

(2.8) (^^)*^^=Idi2(M.). 
Remark. =^^(^^)* ^ Idi2(K2<i). 

Definition 2.2 (Asymptotically high frequency function). Let {it^} C Lp'iW'-) be 
a family of functions such that ||«^||j^2 is uniformly bounded. Given (5 > 0, we say 
that {u^} is asymptotically high frequency with cutoff 5, if 



/ 



\{^'u'){q,p)\' dqdp = 0{en 



lK?'i\Ks 

as £ ^ 0. Ks is the closed set defined in (2.1). 

The definition of asymptotically high frequency functions is motivated by WKB 
function, which is typical in study on high frequency wave propagation. 

Example 2.3. For A{x) G C^{W^) and S{x) e C°^{W^), \VS{x)\ > 5 > 0, the 
family of WKB functions = A(x) exp{iS{x)/e) is asymptotically high frequency 
with cutoff S. 



The Definition 2.2 is also related with the notion of frequency set, microlocal 
support in microlocal analysis. Please refer to [13] for more details. 



CONVERGENCE OF FROZEN GAUSSIAN APPROXIMATION 



5 



3. Hyperbolic system and Hamiltonian flow 

We consider an N x N linear hyperbolic system in d dimensional space, 

d 

(3.1) dtU + Y,Mx)dxiU = 0, 

1=1 

where m = (wi, . . . , wjv)'^ : M'^ and Ai : R'^ ^ M^x^, 1 < I < d sue smooth 

matrix valued functions. 

We assume that the system (3.1) is strictly hyperbolic, i.e., for any {q,p) G 

R^'*, \p\ > 0, the matrix Ylf^i piAi (q) has N distinguished real eigenvalues, denoted 
as {Hn{q-p)}n=i- We denote by Ln{q,p) and Rn{q,p) the corresponding left and 
right eigenvectors, i.e., 

d 

(3.2) ^p;L;^(g,p)^;(g) = H^{q,p)Ll{q,p), 
1=1 

d 

(3.3) J2PlM'l)Rn{q,p) = Hr,{q,p)Rn{q,p), 

1=1 

with the normalization 

Lm{q,p)Rniq,P) = Smn, 

where 6mn is the Kronecker symbol. Note that as Ai{q) is smooth, Hn{q,p), Rn{q,p) 
and Ln{q,p) are smooth functions of {q,p) for |p| > 0. Singularities may occur at 
p = 0. 

The Hamiltonian flow associated with Hn{q,p) is given for |p| > as 
( dQn{t,q,p) 



(3.4) 



= dp„Hn{Qn{t,q,P),Pn{t,q,p)), 

= -dQ^Hn{Qn{t,q,p),Pn{t,q,p)), 



dt 

<iPn{t,q,p) 



with initial conditions 

(3.5) Qn{0,q,p) = q and P„(0,g,p)=p. 

This gives the map (q.p) i-^ {Qn{t,q,p), Pn{t,q,p)) for both t > and < < 
(forward and backward flows). For p = and q £ R'^, we define Qn{t, q,p) = q and 
Pn{t,q,p) =p. 

We make the following assumption for the system (3.1) we consider, which will 
be assumed for the rest of the paper without further indication. 

Assumption A. For each n = 1, • • • ,N, there exists constant C > 0, so that the 
Hamiltonian iJ„ satisfies for any {q,p) G R^*^ with |p| > 

\p-dqHn{q,p)\<C\pf , and \q ■ dpHn{q,p)\ < C \qf . 

Assumption A is understood as a global Lipschitz condition for the ODEs (3.4), 
as we can see from the next Proposition. 
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Proposition 3.1. For each n= 1, - ■ ■ ,N,T>0 and d > 0, there exists a constant 
St > such that {Qn{t,q,p), Pn{t,q,p)) S Ks^, i.e., it satisfies 

(3.6) ST<\Pn{t,q,p)\<l/ST, 

and 

(3.7) \Qn{t,q,p)\<l/ST, 
for any {q,p) & Ks and t € [0,T]. 

Proof. We suppress the subscript n in the proof, since the argument is the same 
for each branch. 

Differentiating \P\ with respect to time and using (3.4) produce 
dt ' ' P dt P ^ 



Hence, by Assumption A, 



^^\P\<T^AP-dQH\<C\P\. 



Gronwall's inequality implies 



m&ji \P{t,q,p)\<CT\p\ 

0<t<T 



On the other hand, differentiating |P| ^ gives 



This yields 



and hence 



max^|P(t,g,p)| '<Ct\p\ V 



Therefore, there exist Ct, Ct > such that for t € [0, T] and any q G M.'^, 
CT\p\<\P{t,q,p)\<CT\p\, 

which implies (3.6). 

The time derivative of \Q\ can be bounded by 

^^\Q\ = ±Q.d,H<C\Q\ 

by Assumption A. Equation (3.7) then follows from Gronwall's inequality. □ 
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Proposition 3.2. For I — 1, ■ ■ ■ ,d, k e M and 6 > 0, there exists constant Ck,s 
such that 

Moreover, for each n = 1, • • • ,N, there exists constant Ck,s such that 

^k,s[Hn] < C'k^S, ^k,s[Rn] < Ck,5, ^k,s[Ln] < Ck,S- 

Proof. The conclusion follows immediately from the smoothness of Ai,Hn, Rn, Ln 
for |p| > and the compactness of the set K^. □ 

We define canonical transformation and action associated with Hamiltonian flow. 



Definition 3.3 (Canonical transformation). Let 
(3.8) K : 



|2d ^ j^2d 



(g,p)^(Q"(9,p),P"(9,p)) 
be differentiable for \p\ > 0. We denote the Jacobian matrix as 



(3.9) i^«(9,p) 



We say k is a canonical transformation if is symplectic for any {q,p) G K^'' with 

bl > 0, 

(3.10) (Fn- ( ' '^A F^^( ' ''A . 

As a corollary to Proposition 3.2, the map Kn.t ■ {q,p) '-^ {Qn{t, q,p), Pn{t, q,p)) 
is a canonical transformation. We formulate this as the next proposition, which 
also gives additional smoothness properties for the Jacobian. 

Proposition 3.4. For each n = I,-- - ,N, the map Kn,t is a canonical transfor- 
mation, and for any k > 0, 6 > and T > 0, there exists constant Ck,5,T such 
that 

(3.11) sup AkAF'"-']<Ck,6,T. 

te[o,T] 

Proof. The argument is the same for different branches, hence we will suppress the 
subscript n for simplicity. 

Proposition 3.1 shows that there exists 6t such that {Q{t, q, p), P{t, q, p)) S Kg^ 
for {q,p) e Kg and t G [0,T]. 

Diff'erentiating F**' {q, p) with respect to time gives 
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Proposition 3.1 and 3.2 imply 

\F^^{q,p)\<C\F^^{q,p)\. 



dpdqH dpdpH ' 
-dgOgH -dqdpH 



where C is independent of {q,p)- Hence, by Gronwall's inequality and F'^'>{q,p) = 
Id2d, one has 

\F-^{q,p)\<exp{C\t\). 
Differentiating (3.12) with {q,p) yields 

The estimate (3.11) follows by an induction argument. □ 

Definition 3.5 (Action). Let k be a canonical transformation defined in (3.8). A 
function S'^ : M^"^ — >■ R is called an action associated to k if it satisfies 



(3.13) dpS^{q,p) = {dj>Q^{q,p))P^{q,p), 

(3.14) d^S^{q,p) = -p+{d,Q^{q,p))P\q,p), 

for any \p\ > 0. 

The action Sn{t,q,p) — S'^"-*{q,p) corresponding to Hn{q,p) solves the following 
equation 

J c 

(3.15) ^ = ^" ■ 9pHn{Pn,Qn) " Hn{Pn,Qn), 

with initial condition Sn{t, q,p) = 0. It is easy to check that 5„ is indeed the action 
associated with Kn,f, which is given by 

(3.16) S„{t,q,p)= I P„{T,q,p) ■ drQ„iT,q,p) - Hn{Q„iT,q,p),P„{T,q,p)) dr. 

Jo 

Now we are ready to construct the Fourier integral operator which will be used 
in the definition of frozen Gaussian approximation in the next section. 



Definition 3.6 (Fourier Integral Operator). For M e L°°(M2d. c^xn-^^ ^ Schwartz- 
class function u G 5(M'^; C^) and n ~ 1, - ■ ■ ,N, we define 

(3.17) {I^{t,M)u){x) = {2ne)-^'^/'^ [ e'*"(*'^'^'«'f)/'^M(g,p)u(y) dgdpdy. 
where the phase function is given by 

(3.18) ^n{t,x,y,q,p) = Sn{t,q,p) + ^ k - Qn{t,q,p)\^ + Pn ■ {x - Qn{t,q,p)) 

+ ^\y-qf -p-{y-q)- 
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Proposition 3.7. If M e L=^{R^'^;C^'''^), for each n = 1,- ■ ■ , N and any t, 
In{t, M) can he extended to a linear hounded operator on L^(IR''; C-'^), and we have 

(3.19) ||T^(t, M)||^(^2(Rd.civ)) < 2 '^^'^ \\M\\j^^^^2d.cNxNy 

Proof. For u,v G L^(]R'^), taking the inner product of v and (3.17) gives 
{v,X^{t,M)u) = (27re)-3'^/2 f ^;(a;)*e**"(*'^'^'«'3')/=M(g,p)u(y)dgdpd2/dx 



= (27r£)-^'*/2 / dqdp, 



,iS„{t,q,p)/e 



:P„-{x-Qr,)-^\x-Qn\ 



v{x) dx 



= 2-<^/2 / e 



xM(q,p) f e-eP<y-''^-hy-i\\{y)dy 
^Sn{t,,,p)/ej^^Q^^ P^)M{q,p){^'u){q,p) dqdp. 



Therefore, 

\{v,I'^{t,M)u)\ < 2-^1^ \\{.9--v) o ^„,J^.(K..^c«) 



L2(R2<i.CW) 

< 2 '^/^ ll^lll,2(]ltd.CN) ||-M"|ll,oo(]R2<i.CiVxN) ||w||j;,2(]R<i.CN) , 

where we have used Proposition 2.1 and the symplecticity of the canonical transform 
K;„^t. This imphes 

ll^n(*)-^)"lll,2(M<i;CJV) < 2 ''^^ ||M|| j^oo(R2<i.cNxiV) j^2(R<i.ciV) , 

which proves the Proposition. □ 
4. Formulation and main result 
Let X5 : '^'^ [0, 1] be a smooth cutoff function for 5 > so that 

1, (9,p) e Ks\ 



0, (g,p) eR2d\^5/. 



and for any k G N , there exists constant Ck.s such that 



sup sup 



< Ck,5- 



The construction of the cutoff function is standard. 

Construct the JC-th order frozen Gaussian approximation with the cutoff 5 as 
(4.1) 



AT 



{n,K,sUo){x) = {K{tMl,K,tX5)Uo) {X) 

1 " [ 

= (27r£)3rf/2 y^3^ e'*"/''A1^_^_4(g,p)x5(g,p)wo(t/)dgdpd?/. 
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Here, the matrix symbol Mn k t given by 

K-l 

(4-2) MIk,Mp) = M^j^it, q,p)=J2 ^''MnAt^ 

fe=0 

Before we introduce the N x N matrix valued function Mn,k{t,Q,p), we need to 
define some operators first. We denote 

(4.3) d^^dq-idp, Zn{t,q,p)=d;,{Qn{t,q,p)+iPn{t,q,p)). 

Note that Zn{t,q,p) is invertible for \p\ > 0, which will be proved in Lemma 5.1. 
Without further indication, we are going to use the Einstein summation convention, 
except for the branch index n. 

We denote b : M?"^ ^ C as a generic function defined on K^"^. For indices 
ii>i2 S {!,••• ,d}, define operators V and Q as 

(4.4) (Vr^jMl,?) = -d.^{h{q,p)Z-]^^{t,q,p)), 

(4.5) iGn,nj2b){q,p) = d^,Q„j^{t,q,p)b{q,p)Z-^.^i{t,q,p), 

provided the right hand sides are well defined (in particular, b is differentiable at 

For z/ e N, € {1, • • • , d} for a = 1, • • • ,u, and £ > 0, define operators T^jj.- 

as 

(4.6) r^p^eVr.,,,b, 

(4-7) rl'J^j^b = egn,nhb + e^Vnj.Vnj.b, 
and for v > 3, 

CL=\ 

provided the right hand sides are well defined at ((7,p) G . Notice that 2?, Q 
and T depend on t as well, which we choose not to make explicit in notation for 
simplicity, with the hope that no confusion will occur. 

It is easy to see that we can rewrite T^lj^...j^ in orders of e as 

(4-9) v:i..,j^= E e^-Ki-.b, 

where [•] is the ceiling function. The last equality defines T^^^...^^- The definition 
of these operators X>, Q and T can be extended to vector-valued and matrix-valued 
functions so that the operators act on the function elementwisely. For example, for 
the matrix-valued function M : M^*^ C^""-^, we have 

{Vn,j,M)kl=Vn,j,MkU 

iov k,l = !,■■■ ,N. 
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We define the operator Cn,k as follows. For M : R x M^'^ £Nxn ^^^^ ^ 

M X R^d ^ith > 0, Cn,k{M) is given by 

(4.10) ) / 

= li^PnjAjiQn) - Hn{Qn,Pn)MN)M{t,q,p), 

{Cn,iM){t,q,p) = dtM{t,q,p) 

+ lDn,j(^{dtPnJ - ldtQn,j)'idN 

(4^1) +iAj{Q^) + P^^idjAi{Q^)jM{t,q,p)^ 

+ lPn,ldljMQn))M{t,q,p)y 

where S„, Qn and Pn on the right hand sides are evaluated at {t, q,p) and we have 
introduced the short-hand, 

9n -j.MQn) = 5q„„vQ„.„^'W-)' any > 1. 
Moreover, for fc > 2, we define 



(4.12) 



2k 

i^n,kM){t,q,p) = Y.V:i-,. {[-T^^d^^y^-AiAQn) 



i'=k 



^Pn,idl...j^Ai{QrS)M{t,q,p) 



We remind the readers that X>, Q and T depend on t implicitly. 

We are ready to give the matrix-valued function M„ ^ now. Mn^k is defined for 
t e [0,T] and {q,p) € R^"^ with \p\ > 0. First, M„,o is given by 

(4.13) M„^o{t,q,p) = an,oit,q,p)RniQn{t,q,p),Pn{t,q,p))Ln{q,p)'^ , 
where (T„,o is determined by the evolution equation, 

(4.14) -^o-„,o(i, q,P) + (^nfi{t, q,p)\n{t, q,p) = 0, 
with the initial condition 

(4.15) a„,o(0,<z,p) = 2-^/2. 
In (4.14), we have used the short- hand 

K{t, q,p) = Ll{dp^Hn ■ dq^Rn - dq^Hn ■ dp^Rn) 

(4.16) - {dM'^{iAj - 5p„,,i?„) + i{dQ„^,Hr. - Pn,idjAi))RnZ-]k 
+ d,,QnjZ-lLl{-djAk + lPn,id%Ai)R^, 
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where Qn,Pn are evaluated at {t,q,p), Aj's are evaluated at Qn, and iJn,L„,i?„ 
are evaluated at {Qn,Pn)- 

Notice that the action £„,o is just multiplication with the matrix t [Pn,j Aj {Qn) — 
Hn{Qn,Pn)^dN) on the left. We define the matrix, 

(4.17) 4,o(9>P) = « E {Hn{q,p) - H^iq,p)y^Rm{q,p)Lliq,p). 
For A; > 1, M„_fc is given by 

(4.18) Mn,k{t,q,p) = an,k{t,q,p)Rn{Qn{t,q,p),Pn{t,q,p))Ln{q,p)'^ 

+ M^,k{t,Q,P), 

where M^^{t, q,p) is given by 

k 

(4.19) M^fc(t, q,p) = Clo{Qn{t, q,p), Pn{t, q,p)) ^(£„,.M„,fe_,)(i, q,p), 

s=l 

and <yn,k{iiQ.iP) solves 

(4.20) + (7„,fcA„ + Ln{Qn{t, q,p), Pn{t, 9,p))^ (^(£„,iM^fe)(t, q,p) 

fe + l N 

+ ^{Cn,sMn,k+l-s){t, q,p) ] Rn{q,p) = 0, 
s=2 ^ 

with the initial condition <Jn,k{0,QjP) = 0. 

We now state the main result. The following theorem indicates that the ii'-th 
order frozen Gaussian approximation (FGA) gives an order 0{e^) approximate 
solution to strictly linear hyperbolic system. 

Theorem 4.1. Consider a strictly linear hyperbolic system (3.1) under Assump- 
tion A. For a family of initial conditions {mq} that is asymptotically high frequency 
with cutoff S, and \\uf^\\ ^2 (jgd^ < M, then for any T > and K £ N, there exist 
constants Ct,k and Eq > 0, such that for e € (0, £o], 

max^ ll^t^o - KK^s^h^r.Hm) < Ct,kMs^ ■ 

Remark. The cutoff 5 is used in the formulation to avoid singularities presented at 
|p| = 0. From an numerical point of view, the cutoff is quite natural. Indeed, in 
the numerical implementation of frozen Gaussian approximation [11,12], a cutoff 
on phase plane is always used for efficiency. 

Remark. The necessity of assumption that the initial value is asymptotically high 
frequency is presented in the following example. It shows that the FGA does not 
work if the FBI transform of the initial condition concentrates around \p\ = as 
£ ^ 0. 
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Example 4.2. Consider the acoustic wave equation in two dimension with constant 
coefRcients, 

dfU + AidxiU + A2dx^u = 0, 

where 





/o 









/o 





























V 











1 





Then the eigenfunctions in (3.2)-(3.3) arc given by 

H,{q,p)=0, H±iq,p) = ±\p\, 

where we have used the subscripts ± instead of number subscripts. This imphes 
the system is strictly hyperbohc, and the corresponding eigenvectors are 




R±{q,p) = I ±P2 
\P\ 



p2 \ ^ I ±pi/ br 

Liiq,p) = -"72 I -Pi , L±{q,p) = - ±^2/ \p\" 



\P\ 







1/W 



It is easy to see that the system satisfies Assumption A. 

We focus on the "+" branch and omit the subscript "+" for convenience, then 
the Hamiltonian flow is given by 

r dg ^ p_ 

dt \p\' 

dP 



(4.21) 



dt 



0, 



with the initial conditions Q(0,q,p) = q and P{0,q,p) = p, which implies the 
following solution, 



(4.22) 

Therefore 

(4.23) 



P 



Q = q + f-t, P = p. 
\P\ 



detZ = 4-^. 

\P\ 



According to the equation (2.7) in [12], the solution of amplitude a{t,q,p) is 



(4.24) 



2it 

¥\ 



1/2 



with the branch of square root determined continuously in time by initial condition. 
If we take the initial condition of u as 



(4.25) 



uo{x) = (uo,i,Mo,2,uo,3)'^ = (exp(-|a;|^/2), 0, 0)^ 



14 



JIANFENG LU AND XU YANG 



then the FBI transform gives 



■ exp 



(4.26) 



ip-q 



\P\' 



1 + s 2(1 + e) 2s{l + e)J' 



(7r£)V2(l+£) 
^'(W0,2) = 0, 

The leading order frozen Gaussian approximation of u is given by 
(4.27) u^t,x) = ^{^r{cT{t,q,p)R{Q,P)L''{q,p)^'{uo)), 



By Proposition 2.1 
(4.28) 

1 



a{t,q,p)R{Q,P)L'{q,p)^'{uo) 



L2 



Notice that, as £ ^ 0, 



it 



A. 

4\p\< 



exp 



\P\' 



1 + e e{l + e) 



dq dp. 



I 



1^1 



e-T+7d9 = C»(l), 



and 

1 







/ 


2bl 



pi _ , 1 

e ^rr+;) dp = - 



A\p\ 



OO /•27T 

Jo 



2r 



cos 



-e r dr d0, 





^ it 






Jo Jo 




-1/2)^ 





cos'' 



e TT+JTdrd6' 



where the second equahty is obtained by change of variable. 



Hence one has 



lL2 



is of the order 0{£ ^/''), while 



lL2 



under the initial 



condition (4.25) is of the order 0{1), which implies FGA can not be a good ap- 
proximation to the acoustic wave equation under this choice of initial condition. 



5. High order approximation 

In this section we introduce a high order approximation to the solution of (3.1) 
based on frozen Gaussian approximation. This is a key step in the proof of Theo- 
rem 4.1. We state and prove some preliminary lemmas first. 

For a canonical transformation Kn,t, define Z'^'^'*{q,p) for \p\ > as 

{q, p) = (0'=-' (g, p) + zP«".' (g, p)) , 
which is related to Zn{t,q,p) defined in (4.3) by 



Z„(t,g,p) = Z«-*(g,p). 
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Lemma 5.1. Z'^'^'*{q,p) is invertible for each n = 1, • • • ,N and {q,p) G M^'^ with 
\p\ > 0. Moreover, for any k>0 and S > 0, there exists constant Ck,s such that 

(5.1) Afc,5[(Z«".*)-i(g,p)] <C7fe,5. 

Proof. Since the proof is the same for each branch, we omit the subscript n in 
notation for convenience. 

Observe that Z'^*{q,p) can be rewritten as 



Z^'{q,p) = d4Q^'{q,p)+iP'^^{q,p)) = {tldd Id^) {F'^^f{q,p) 
Therefore, 



-ildd 
Idd , 



= (tldd ldd){iF^^fF^^)iq,p) 



-ildd 



= [zldd Idd) {{F^^fF^^){q,p) (^-^^^'^ + 2ldd. 

In the last equality, we have used the property of symplecticity of F'^*{q,p). There- 
fore is positive definite, which implies det(Z'^*((7,p)) is bounded 
away uniformly from zero for \p\ > 0. The invertibility of and bounds (5.1) 
follows from the representation of {Z'^*)~^{q,p) by minors and Proposition 3.4. □ 

The following Lemma plays an important role in frozen Gaussian approximation. 

Lemma 5.2. For each n = I, ■ ■ ■ , N , t € [0,T], let b{y,q,p) : R^d _> such that 
for any y € R'', supp6(y, •, •) C Kg. Assume that there exists some m e M that for 
any fc G N, 

sup(l + |2/|')-"'/2Afc,,-[6(2/,-,-)] <oo. 
Then for a multi-index (j'l, • • • , jV) with u > 1, 

(5.2) (a; - Qn{t,q,p))h ■■■{x- Qn{t,q,p))j^h{y,q,p) ~ T^^-^...j^b{y,q,p). 

Here T^'j^...j^ is defined by the recursive relation (4.6)-(4.8), corresponding to the 
n-th branch. Here, we have used the notation f g, if 

(5.3) / f{y,q,p)ei^-^''-'y''''P^dydpdq= [ 5(2/, g,7^)e**"^*'"'"'*'^^ dy dpdg. 
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Proof. We omit the subscript n in the proof for simphcity, because the argument 
is the same for each branch. 

Observe that at t = 0, for \p\ > 0, 

dgS{0,q,p) - dgQ{0,q,p)PiO,q,p)+p = 0, 
dj,SiO,q,p) - dpQ{0,q,p)P{0,q,p) = 0. 

Using (3.4) and (3.15), we have 

-^(9,5 - dgQP + p) = dgidcS) - dg{dtQ)P - dgQdcP 

= dq{P- dpH -H)- dq{dpH)P + dqQdqH 

= dgPdpH - {dgQdq + dqPdp)H + dqQdQH = 0, 

where 5, Q and P arc evaluated at (t, q,p), and dqH, dpH are evaluated at {Q, P). 
Analogously we have 

^ {dj,S{t, q,p) - dpQ{t, q,p)P{t, q,p)) = 0. 

Therefore for all t G [Q,T], we have 

(5.4) dgS{t, q,p) - dgQ{t, q,p)P{t, q,p)+p = 0, 

(5.5) dpS{t, q, p) - dpQ{t, q, p)P{t, q, p) = 0. 

Straightforward calculations yield 

dq^{t,x,y,q,p) = {dqP{t,q,p) - idqQ{t,q,p)){x - Q{t,q,p)) -i{y-q), 
dp^{t,x,y,q,p) = {dpP{t,q,p) - idpQ{t,q,p)){x - Q{t,q,p)) - (y-q). 

This implies that 

(5.6) id^^{t,x,y,q,p) = Z{t,q,p){x - Q(t,q,p)), 

where dz and Z{t, q,p) are defined in (4.3). 

Using (5.6) and Lemma 5.1 for the invertibility of Z{t,q,p) on the support of 
b{y^ 9,P), one has, for j = 1, • • • , d, 

[ {x- 0(i,(j,p)) et*(*'-'^'«.f)6(y,g,p)dz/dgdp 

= e f ZT^{t,q,p) (-dz,^{t,x,y,q,p)) et*(*'-'^'«'f)6(y, g,p) dy d^dp. 

We remark that the integrability of the above integral follows from the exponen- 
tial decay in y of ex.p{i^{t,x,y,q,p)/e) and compact support of b{y,q,p) in {q,p). 
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Integration by parts gives 

f {x- Q{t,q,p)).eT^^''^'y''i'PMy,q,p)dydqdp 

= -e [ d,,{Z-\t,q,p)b{y,q,p))ei^^'''''y''i'P^dydqdp 

= [ £Dj6(y,9,p)et*(*'^'2''«'f)dydgdp 

= / 7;.''"6(y,g,p)et*(*'-'^'«'f)dydgdp. 

This proves (5.2) for v = 1. 

Making use of the above equaUty twice produces (5.2) for u = 2, 

{x - Q{t,q,p)).^ {x- Q{t,q,p))^^b{y,q,p) 

~ eVjJ{x - Q{t, q,p))_.b{y, q,p) 



= £Qjij2Ky,q^P)+s{x - Q{t,q,p)) _.Vj^b{y,q,p) 
= Tj;lb{y,q,p). 



Furthermore, we have 



(a; - Q{t,q,p)).^ ■ ■ ■ {x - Q{t,q,p)) .^b{y,q,p) 



'31 

- Q(t,<l^P))j^ • • • (x - Q{t,q,p)) _^^ b{y,q,p)^ 
e{x-Q{t,q,p)).^ ■■■{x- Q{t,q,p)) .^_Vj^b{y,q,p) 



+ e^{x- Q{t,q,p))^^ ■■■{x- Q{t,q,p))_.^_^ 

a=l 

X {x- Q{t,q,p)).^^^ ■■■{x- Q{t,q,p)).^_^gj^jJ{y,q,p) 

T-;l.jAv,(i,p), 



where the last step follows by the recursive relation (4.8). This proves (5.2) for 
v>2>. 

□ 



To construct a high order approximation to the solution, we need a filtered 
version of frozen Gaussian approximation, defined as follows. The relationship 
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between ^ g and ^ s ^^^^ ^® addressed in Lemma 5.6 in tlic end of this section. 
(5.7) 

N 

^ ^ r ~ 
= (2^^)3d/2 E y^^/'*''^*'^'^''''^^/^A^^,K,M(9,p)«o(2/)dgdpdy. 

Here the matrix- valued function A4^ Kts- ^^"^ ~^ (^nxn jg gjygj^ 

(5.8) M'^^K,tAl^P) = M'^^kA*' 1^p)=Y1 ^^Mn,kAt^ <1'P)- 

fe=0 

For {q,p) € IR2<^ and t € [0,r], M„^k,sit,q,p) is defined as follows. For k = 0, 

(5.9) Mn.o,s{t, q,p) = ^nfl,s{t, q,p)Rn{Qn(t, q,p),Pn{t, q, p)) L]^ (q, p)xs{q, p) , 
where xsiltP) is the filter function given before, and a^n,o,s{t,q,p) satisfies 

(5.10) -^'^n,o,5{t,q,p)X6{q,p)+'^nfi,5{t,q,p)\n{t,q,p)X6{q,p) = 0, 

with the initial condition ct„ n ,s{0,q,p) = 2'^/^ and Xn{t,q,p) is given by (4.16). 
For fc > 1, Mn^k,S is given recursively by 

(5.11) Mnxs{t,q,p) = an,k,5{t,q,p)Rn{Qn,Pn)Ln{q,p)'^Xs{q,p) + M^,^s{t,q,p), 
where M^,. g{t,q,p) is given by 

fc 

(5.12) M^fc_5(t,g,p) =4^o(Q„,P„)^(£„,,M„,fc_,,5)(i,g,p), 

s=l 

and <Jn,k,s{t,q,p) satisfies 

-^'^n,k,5(t, q,p)xs{q,p) + 5„,fc,5(t, q,p)Xn{t, q,p)xs{q,p) 

(5.13) + LniQn, PnV (^{jCn,lM^,k,s)i*' I^P) 

fe+1 N 

+ ^(>C„,sM„,fe+i-s,5)(t, q,p) j Rn{q,p) = 0, 

with the initial condition ^n,k,s{0,q,p) = 0. In (5.12) and (5.13), Qn and P„ are 
evaluated at {t,q,p). We remark that, from the definitions (5.9) and (5.11), it is 
clear that the vahic of an^k,s{t,q,p) outside suppxs{q,p) C will not affect the 
value of MnX6{t,q,p)- 

Lemma 5.3. For each n = 1, - ■ ■ ,N, k G N and any S > 0, there exists constant 
Ck,s that 

(5.14) Afc,44^o(g,p)] < Ck,s- 
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Proof. Recall that £„ o{q,p) is defined in (4.17) as 



'^ifiil^P) = 1 XI {^n{q,P) - Hm{q,p)) ^Rm.{q,P)L]n{q,p). 



By strict hyperbolicity and compactness, there exists gs > such that 



(5.15) 



inf inf \Hn{q,p) - Hm{q,p)\ > gs- 



The estimate (5.14) follows easily. 



□ 



Lemma 5.4. For each n = I,-- - ,N, any t £ [0,T] and any k G N, we have 
suppMn,k,5it,-,-) C Ks/2, suppdtMnxsit,-,-) C Ks/2- Morcovcr, for any s e N, 
there exists constant Ck,s,5,T that 



Proof. The argument is the same for each branch, so we will omit the subscript n 
in the proof. 

By Proposition 3.1, there exists constant 5t such that for t G [0,T], we have 
{Q{t,q,p),P{t,q,p)) e for {q,p) G Ks/2- Prom Lemma 5.1 and Proposition 3.2, 
it is easy to conclude that for any s e N and I = 0,1, 

(5.16) sup A,^s/2[dlX{t,q,p)]<Cs,s,T- 

te[o,T] 

Equation (5.9) implies that 

suppMo,5(i,9,p) C suppx5(g,p) C Ks/2- 
Moreover, (5.10) implies that 

-^\'^o,s{t,q,p)x5{q,p)\ < \Kt^Q^P)\\^o,6{t,q,p)xs{q,p)\ - 
Hence Gronwall's inequality yields 



Taking derivatives of (5.10) with respect to {t,q,p), we have for multi-index a and 



sup f^s,5/2Wn,kAt,q,P)] < Ck,s,S,T\ 

te[o,T] 



sup A^_5/2[5t-W„,fe,d-(f,g,p)] < Ck,s,s,T- 
te[o,T] 



sup sup \(Jo,5{t,q,p)X5{q,p)\<C5^T- 

te[o,T] (q.p)eKs,2 



Z = 0,1, 
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Using (5.16) and by induction, wc have for ? = 0, 1, 



sup K,S/2[dl{cro,5{t,q,p)x5{Q,P))] <Cs,s,t, 
te[o,T] 



which implies 



(5.17) sup A,^s/2[dlMo,6{t,'l,P)]<Ce,S,T. 

te[o,T] 

Next, we consider Mi^s{t,q,p). (5.12) gives 

M^^sit, Q,P) = J^l{Qn{t, q,p), Pn{t, q,p)) {CiMo,s){t, q,p). 

As snpp Mo^s{t,q,p) C Kg/2, we have supp M^g{t,q,p) C Ks/2- By Lemma 5.3, 
(4.11) and (5.17), it is not hard to see that for I = 0,1, 

sup As.s/2[dlM^^s{t,Q,P)] < Cs,s,T- 

te[o,T] 

Similarly, we get from (5.13) the estimate 

sup K,5/2[d\'^iAt^^^P)X&{q,P)] < Cs,5,T- 
te[o,T] 

In summary, we obtain suppMi^5(i, q,p) C Kg/2: and for I = 0,1, 

sup As,s/2[dlMiAt,q,p)] < Cs,5,t- 

te[o/r] 

Continuing this procedure, we obtain the estimates for higher order asymptotic 
terms: supp Mk,s{t,q,p) C supp 7^5/2 for k>0, and 

sup A^^s/2[Mk,s{t,q,p)] < Ck,s,S,T, 

te[o,T] 

sup A^^s/2[dtMk,s{t,q,p)] < Ck,s,s,T- 

te[o,T] 

This proves the Lemma. 

□ 

We now show that the filtered frozen Gaussian approximation gives a high order 
approximation to the solution. 

Proposition 5.5. Under the same assumption of Theorem 4-1, for any T > 0, 

K gN, there exists constant Ct,k, so that for any £> 0, 



{dt + Ai{x)d^;)VlK,sul 



K-1 



Proof. Without further indication, Sn, Qn and P„ are all evaluated at {t,q,p) in 
the proof. 
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For each I = 1, - ■ ■ ,d, Taylor expansion of Ai{x) around x — Qn gives 
2K ^ 

(5.18) Ai{x) - Ai{Qn) + ^ ^(a^ - Qn)n ' ' ' (^^ - Qn)j^d^,...j^Ai{Qn) 



(2_ftr)! ^"^ ■ ■ ■ (-^ Qn)j2K+i-^l,ji---j2K+i{Qn)j 



where ^iji...j2i(:+i(<3n) is the integral 

(5.19) Ai,,,...j,^^,{Q^) = f\l - Tf''d]^+l,^^Ai{Qn+T{x - g„)) dr. 



We define the operators £„,fc on Ci([0, T], C°°(R2<i; C^xJV)) follows. For M : 
we define iov < k < K, 



{Cn,kM){t,q,p) = {jCn,kM){t,q,p), 
iov K + 1< k <2K + 2, 



~ /I 
(£„,.M)(i,g,p) = ^r„-t..,J (-7;;3T)T^r4-.%(Q») 

^2/f+l,ft f ( 1 Q2if ^ 

^ 'n,ji---j2ft:+i I ^2/^)! ■'^■■■J^J*^ ■'^^+^^ ' 



■■j2K+l 

{Qn))M{t,q,p) 



{2K)\ 

+ '^^,j^--j2K+2 (^'(^^]^j2K + 2,3l-j2K + l {Qn)M{t, q,p)j , 

and for k>2K + 3, 

{Cn,kM){t,q,p)=0. 

Differentiating ^n{t,x,y,q,p) with respect to t and x gives 

at$„(i, a;, J/, q,p) = dtS^ - Pn ■ dtQn + {x - Qn) ■ {dtPn - ^^tQn), 
dx^n{t,x,y,q,p) =l{x-Qn)+Pn- 

Combining these expressions, substituting (5.7) and (5.18) into (3.1), and orga- 
nizing terms in orders of s produce, after straightforward calculations, 

fe=0 
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where we have 

N 



Vt,iA^) = $^(27r£)-3'^/2 / ((£„,oM„,i,5)(t, q,p) 

+ (£„,iM„,o,5)(t,g,p))e'*"/^ug(y)dydpdg, 

vt,k,5{x) - V(27r£)-3'^/2 / V(£„,«M„,fe_«,5)(t,g,p)e'*''/^ug(2/)d2/dpd9, 
and the remainder, 

2K+1 K-1 N 

n,KA^) = E E E(2-^)"''^' 

fe=0 s=0 n=l 

X / (£n,fe+s+l M„,;f_s_i,5)(i,g,p)e'*"/^ug(?/) dj/dpd?. 

Here is evaluated at {t,x,y,q,p). 

We will show that Vt,k,S = for fc = 0, • • • , — 1 and control the norm of the 
remainder rt,K,s to prove the proposition. 

By definitions of Rn{Qn,Pn) as in (3.3) and of 5„ as in (3.15), we get 

(dtSn - Pn ■ dtQn + Pn,lMQn)) RniQn, Pn) 

= (dtSn - Pn ■ dtQn + Hn{Qn, Pn)) Rn{Qn, Pn) = 0, 

which implies (£„^o-^n,o,(5)(i, 9jP) = for each n, since M„_o,5 is given as (5.9). 
Therefore, vt.0,5 = 0. 

To prove vt,i,s = 0, it suffices to show that, for each m = 1, • • • , A'', 



(5.20) Ll^iQn, Pn) (^(£n,oM„,i,5)(t, q,p) + (£„,iM„,o,5)(t, q,p)j = 0. 
When n, (5.20) is equivalent to, by (3.3), 

(5.21) Ll{Qn,Pn)({CnfiM^,l,s)i*'1^P) + (£n,lM„,0,5) (t, 5, p)) = 0, 

which is valid by (4.17) and (5.12). 

When m = n, (5.20) is equivalent to, by (3.2), 

(5.22) Ll{Q„, P„)(£„,iM„,o,5)(t, q,p) = 0, 

which will be proved as follows. For convenience, we will omit the index n in the 
following calculations. Substituting (5.9) in (5.22) and using (4.4), (4.5) and (4.11), 
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we can rewrite (5.22) as 

L'^iQ, P) (^dt{ao,s{t, q,p)RiQ, P))L'^{q,p)Xs{q,p) 

- d,, ({idtPj + dtQj - Aj{Q) + iPidjMQ)) 

(5.23) X aoAt,q,P)RiQ,P)Z-kL^iq,p)X6{q,p) 
+ d^^Qji-djAkiQ) + ^Pid^MQ)) 

X ao,s{t,q,p)R{Q,P)Z^J^L'^{q,p)xs{q,p) ) = 0. 
For the first term, easy calculations yield 

L'^iQ, P) (dt{ao,s{t, q,p)R{Q, P))L^{q,p)xs{q,p)] 

[dtaoAt, q,P) + ^oAt, q,P)L'^{Q, P) (dpH{Q, P) ■ dQR{Q, P) 

- dQH{Q,P) ■ dpR{Q,P))^L'^{q,p)xs{q,p). 

To simplify the second term in (5.23), we observe that differentiating (3.3) with 
respect to P and Q yields, 

MQ)R{Q,P) + P,AjiQ)dp,R{Q,P) 

= dp,H{Q, P)R{Q, P) + H{Q, P)dp,R{Q, P), 

and 

PjdiAj{Q)R{Q, P) + PjAj{Q)dQ,R{Q, P) 

= dQ,H{Q, P)R{Q, P) + H{Q, P)dQ,R{Q, P). 
Taking inner product with L{Q, P) on the left produces 

(5.24) L'^iQ, P){Ai{Q) - dp,H{Q, P))R{Q, P) = 0, 

(5.25) L'^{Q,P){PjdiAj{Q) - dQ,H{Q,P))R{Q,P) = 0. 

Define the short hand notation 
Fj{t, q,p) = [idtPj + dtQj)R{Q, P) - Aj{Q)R{Q, P) + tPidjAi{Q)R{Q, P) 

= -(^{A,{Q) - dpM{Q,P))+i{dQM{Q,P) - PidjAi{Q)))R{Q,P). 

Using (5.24) and (5.25), it is clear that for any j 1, . . . , d, 

L^{Q,P)Fj{t,q,p) = Q. 

Hence, 

L^{Q,P)d,,{{ao,sFjZj^){t,q,p)L^{q,p)X&{q,p)) 

= aoAt,q,P)L'^{Q,P) [dz, Fj ZJ^ ){t,q,p)L^{q,p)x5{q,p) 

= -aoAt,q,P){dz,Lf{Q,P) {Fj ZJ^ ){t,q,p)L^{q, p)xs {q,p). 
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Therefore, (5.23) is implied by 

dt^o,s{t,q,p)xs{q,p) + ^o.5it,q,p)L^{dpH ■ dqR - dqH ■ dpR)xs{q,p) 

+ ao,s{t,q,p) (d,^ L)'^{Fj ZJ^) {t,q,p)xs{q, p) 

+ 5o,5 (i, q, p)dz, Qj Z^^ {t,q,p) 

X L^{-d,AkiQ) + '-Pid%Ai{Q))xs{q,p) = 0, 

which is vahd by (5.10) and (4.16). Here in the last equation, H, L and R are 
evaluated at {Q,P). 

Similarly, we can show that Vt,k,s{x) = for higher orders, by using (5.11)-(5.13). 

Therefore, to prove the Lemma, it remains to bound the remainder rt,K,s{x)- 
Using Proposition 3.7, it suffices to control the L°° norm of {Ln,kMn,s,5){t, q,p) for 
k = 1, - ■ ■ , 2K + 2 and s = 0, • • • ,K — 1. This in turn follows from Lemma 5.4 and 
the definition of Ln^k- O 

To relate 'P^ g with g defined in (4.1), the following lemma shows that they 
are essentially the same when applying on asymptotically high frequency initial 
data. 



Lemma 5.6. For any T > Q, K eN, 

sup VIk,su'o - 'PIk,sUo 

0<t<T 



I,2(E<i;CiV) 



Proof. By definition (4.1) and (5.7), 

^ ' n=l k=0 

X s'' [Mn,k{t,q,p)xsiq,p) - Mn,k,s{t,q,p)^uo{y) Aqdpdy. 

From the constructions of M„_fe and Mn^k,6, it is easy to see that for t G [0, T] and 
{q,p) e Kg. 

(5.26) Mn,k{t, q,p) = Mn,kAt, Q,P)- 
As Xs{q,p) = 1 for {q,p) € Ks,we have 

Mn,k{t,q,p)x5iq,p) = Mn,k,s{t,q,p)- 
Using (5.26) with 6/2, we have then 

(5.27) Mn,k{t,q,p) = M^,k,s/2it,q,p) 

for {q,p) G Ks/2, and hence in particular, for {q,p) S suppxa- Combining (5.27) 
with Lemma 5.4 gives 

(5.28) sup sup \Mn,k{t,q,p)xs{q,p)\ < Ct- 

te[0,T] (g,p)eR=<* 
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Lemma 5.4 guarantees supp Mn,k,s{t, ■, •) C Kgf2- Since suppx^ C Kg^2, we also 
have svLpp{Mn,k{t, ■ r)xs) C Ks/2- Therefore, by (5.26), we have 



(5.29) supp(^M„,fe(t, •, ■)xs - M„^k,s{t, •, •) j C Ks/2\Ks. 

Using a similar argiimcnt as in the proof of Proposition 3.7, one has 



L2(R2d.CN) 



N K~l 

< 2"'^/' EE II {MnAt, ■)XS - Mn,kAt, 

n=l fe=0 

TV if-l 

< 2"'/' E E •' •)^^ - Mn,k,5{t, ; •)) 

n=l fe=0 

X ||=^^Wollz,2(Kj/2\^fii;CJ^) 

< C'5,T,if l|-^''WollL2(K5/2\ifj;CN) ' 

where we have used (5.29) in the second inequality, and (5.28) and Lemma 5.4 in 
the last inequality. The proof is concluded by noticing that 

by Definition 2.2. □ 
6. Proof of the main results 

We recall the energy estimate for linear strictly hyperbolic system. The proof 
can be found for example in [21]. 

Lemma 6.1. Given strictly hyperbolic system 

d 

dtu + ^Ai{x)dx,u = /, 
1=1 

with initial condition u{0,x) = uo{x), where Ai{x) are given as in (3.1). For any 
T > 0, there exists a constant Ct such that 

i-T 



sup ||u(t, x)||^2(Rd c«) < Ct( ||Mo(a;)||i2(Rd.civ) + / ||/(s, 

0<t<T V JO 



^)llt2(Rd.CN) ds I . 



Proposition 6.2. Under the same assumption of Theorem 4-1, for any K G N, 
T > 0, there exists constants Ct,k and eq > 0, such that for any e e (0, sq], 



(6.1) 

Proof. We denote 



sup 
te[o,T] 



e'{t,x) = {Vtul){x) - {VlK,snl){x). 



Proposition 5.5 implies 



sup \\{dt+ Ai{x)d^i)e^{t,-)\\j^.2,^i.cN^ < CT,KMe 
te[o,T] 



K-l 
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Notice that by the construction of filtered frozen Gaussian approximation, 

X {q, p)Ll {q, p)xs {q, pXo (y) dq dp dy 

where we have used that fact that for {q,p) S M^'' with |p| > 0, 

X] ^n{q,p)Ln{q,p) = Idjv. 

n=l 

This implies 

e/(0,x) - ulix) - {K,K,6ul){x) = - Xs)J''ul). 

Hence, using Proposition 2.1, 

Ik^CO, Olli^^Rd-ciV) = 11(1 - X(5)-^^WollL2(M2<i;C«) ^ \\^^'^l\\L'^{9?d\Ks;C'^) = 0{e°°). 

The conclusion of the Proposition follows easily from Lemma 6.1. □ 

Finally, we conclude with the proof of Theorem 4.1. 
Proof of Theorem 4-1- Triangle inequality gives 



+ 



L2(E<i;C«) 



+ \\'^t,K+l,S'^0 -'^tV,5"o|li,2(K<i.cN) • 

The first two terms are estimated by Lemma 5.6 and Proposition 6.2. For the last 
term, notice that by definition 

Kk+^suI - VIk,5uI = I^(i,£-^M„,K(t, •, ■)x5Xo, 
and hence, using (5.28) and Proposition 3.7, we have 

ll^tV+l.'S^O - ^tV,5"o|li2(R<i.ciV) < Ck,tM£^. 
The Theorem is proved. □ 
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